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ABSTRACT 



Context. The blazar sequence is a scenario in which the bolometric luminosity of the blazar governs the appearance of its spectral 
energy distribution. The most prominent result is the significant negative correlation between the synchrotron peak frequencies and 
the synchrotron peak luminosities of the blazar population. 

Aims. Observational studies of the blazar sequence have, in general, neglected the elfect of Doppler boosting. We study the depen- 
dence of both the synchrotron peak frequency and luminosity with Doppler-corrected quantities. 

Methods. We determine the spectral energy distributions of 135 radio-bright AGN and find the best-fit parabolic function for the 
distribution to quantify their synchrotron emission. The corresponding measurements of synchrotron peak luminosities and frequen- 
cies are Doppler-corrected with a new set of Doppler factors calculated from variability data. The relevant correlations for the blazar 
sequence are determined for these intrinsic quantities. 

Results. The Doppler factor depends strongly on the synchrotron peak frequency, the lower energy sources being more boosted. 
Applying the Doppler correction to the peak frequencies and luminosities annuls the negative correlation between the two quantities, 
which becomes positive. For BL Lacertae objects, the positive correlation is particularly strong. 

Conclusions. The blazar sequence, when defined as the anticorrelation between the peak frequency and luminosity of the synchrotron 
component of the spectral energy distribution, disappears when the intrinsic, Doppler-corrected values are used. It is an observational 
phenomenon created by variable Doppler boosting across the synchrotron peak frequency range. 

Key words, galaxies: active - BL Lacertae objects: general - quasars: general - radiation mechanisms: non-thermal 



1. Introduction 

The most powerful active galactic nuclei (AGN) are active emit- 
ters at almost all frequencies and exhibit both strong flaring and 
typically high polarization. These sources are referred to, in gen- 
eral, as blazars. The source of their power is a supermassive 
black hole embedded in the host galaxy, which generates a mas- 
sive accretion process. The bulk of the continuum radiation from 
these sources does not, however, originate in the nucleus, but 
in two jets that emerge symmetrically from the core. A tell-tale 
sign of the presence of jets is strong radio continuum, which is 
believed to be produced via synchrotron process by relativistic 
electrons spiralling in the magnetic field of the jet. This syn- 
chrotron radiation can be seen in the blazar spectral energy dis- 
tribution (SED) as a bump at radio to X-ray frequencies in the 
log V - log (vFy) -representation. It is followed by another bump, 
which is often attributed to the inverse Compton process. 

The shape and position of the synchrotron SED component 
in the log (vFy) versus log v coordinate system and the mech- 
anisms that control its properties have been discussed for a 
number of years. It was first suggested by Padovani & Giommi 
(i 1995b that the peak of the synchrotron component can occur at 
any frequency from radio to soft X-rays, producing detectable 
differences in the source properties, such as broad band spec- 
tral indices. The position of the synchrotron component peak 
on the log v -axis, the synchrotron peak frequency Vp, was 
linked to the inttinsic s ynchrotron peak lumin osity, VpLy^p, by 
iFossati et alj (Il998l) and lGhiselUni et alj (Il998h . They proposed 



that there was a negative correlation between the two quanti- 
ties, that is sources with their peak frequency in the radio band 
would be more luminous than those that had their peak in the 
X-ray band. This correlation would be created by decreasing in- 
trinsic power and decreasing cooling effects, caused by exter- 
nal radiation field, as the sync hrotron peak moves from low to 
high energies (Ghiselli ni et aIJ [T998). The blazar luminosity se- 
quence has been studied by many authors (e.g..lPadovani et akl 
120031: ICaccianiga & Marchall2004l: [Anton & Brownell2005L for 
a concise review of the present status of the blazar sequence, see 
Padovani ( |2007|) ). Lately, evidence against it has been mount- 
ing; as the number of blazars known has increased, it has been 
found that the previo usly-observed coiTelation does n ot describe 
the newer data sets. iGhisellini & Tavecchiol (l2008h attempted 
to explain these outliers using a revised theoretical blazar se- 
quence, which related the SED shape to the mass of the central 
black hole and its accretion rate. This new scheme reduced the 
tightness of the correlation between the peak frequency and lu- 
minosity, allowing the presence of high-luminosity high-Vy, ob- 
jects, and even predicting the existence of low-power low-Vp 
objects. However, the negative trend be tween the two quanti- 
ties was, in general, expected to hold. In lNieppoIa et alj ( l2006l 
hereafter Paper I) the synchrotron peak frequencies of more than 
300 BL Lacertae objects (BL Lacs, BLOs) were determined and 
used to plot the VpLy p versus Vp -correlation. No significant anti- 
coiTelations could be found using this extensive sample, which 
included many high-energy sources. 
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One important matter overlooked in these observational 
studies of the blazar sequence was the effect of Doppler 
boosting. When the source of the continuum radiation (i.e., 
the shocked plasma in the jet) is moving towards us, the 
observed radiation is blue shifted and enhanced significantly 
jBlandford & Konigilll979l) . The amount of Doppler boosting 
depends on the viewing angle of the jet, 9, and the Lorentz 
factor y. It is therefore clear that the Doppler factor D is 
by no means similar for all AGN, or even for all blazars, 
which implies that it can significantly eff'ect the appearance 
of the blazar sequence. In the original bl azar sequence papers 
jFossati et al.' 1998; GhiselHni et al. 1998) and in theoretical 
work (Ghisellini & Tavecchio 2008), the Doppler factor was as- 
sumed to be constant irrespective of the peak frequency Vp. In 
this paper, we demonstrate that this assumption is inaccurate. We 
examine the Vp versus VpLy p -correlation for a complete sample 
of AGN using Doppler-corrected values based on new estimates 
of Doppler-boosting factors (Hovatta et al. 2008, in preparation). 
In this way, we are able to determine if there is any pattern in the 
intrinsic, true properties of the sources. 

The paper is organized as follows: in Sect. |2] we present the 
source data that we study, and in Sect.|3]we describe the methods 
of acquiring the data. In Sect. |4] we present the results, and fin- 
ish with a discussion and conclusions in Sects. |5] and |6] respec- 
tively. Henceforth, all luminosities are frequency dependent. For 
clarity, we denote luminosities simply with L instead of Ly, and 
use the subscript p to denote peak frequency and luminosity of 
the synchrotron component determined directly from the SED. 
Throughout this paper, we assume Hq = 65kms"' Mpc"' and 

2. Sample 

We study sources selected from the Metsahovi Radio 
Observatory source Ust, which meet one or both of the following 
criteria: 

1 . The source has a declination of > deg and the highest ob- 
served flux exceeds 1.2 Jy at 37 GHz. 

2. The source has a Doppler fa ctor determined by either 
iLahteenmaki & Valtaoial d 19991) or Hovatta et al. 2008 (in 
preparation). 

Sources meeting the first criterion comprise a complete 
northern flux-limited radio-loud AGN sample; they constitute 
a large majority of the entire sample, 113 sources out of 135. 
The remaining sources are radio loud as well, but do not reach 
the flux limit of 1.2 Jy, or are southern sources. They are in- 
cluded because they have well sampled radio flares at 22 or 
37 GHz, allowing a reliable determination of the Doppler fac- 
tor from variability data. The full sample consists of 34 highly 
polarized quasars (HPQ), 33 quasars of low polarization (LPQ), 
31 BL Lacertae objects (BLO), 26 quasars with no polarization 
data (QSO), 9 radio galaxies (GAL), and 2 unclassified objects. 
HPQs have a measured polarization of > 3 % in the literature 
and LPQs have a polarzation of < 3 %. BLOs are included in 
IVeron-Cettv & Veron (2000), or defined as such in the litera- 
ture. GALs are radio galaxies of varying properties, which all 
are non-quasars. 

The full sample is listed in Table[T] Columns (1) and (2) pro- 
vide alternative names for the source and Cols. (3) and (4) in- 
dicate the right ascension and declination, respectively. Column 
(5) indicates whether the source belongs to the 1.2 Jy complete 
sample. Column (6) provides the classification of the source, our 



measeurement of the synchrotron peak frequency (see Sect. 13. It 
is listed in Col. (7), and Col. (8) gives the redshift of the source. 
Finally, Col. (9) indicates the sources that belong to the blazar 
sample discussed in Sect. 14.21 

3. Methods 

3.1. Spectral energy distributions 

The spectral energy distributions of our selected sample were 
determined in a similar manner to those of Paper I. A sig- 
nificant quantity of multi-frequency data were collected from 
various sources, including th e CATS databas jj of th e Special 
Astrophysical Observatory (IVerkhodanov et alTl2005b and the 
NED databas^ at IPAC. Most sources have a low-energy syn- 
chrotron peak and we were able to find a substantial amount of 
radio data. The majority of the high frequency ra dio data were 
acqui red from Metsahovi Radio Observatory (ISalonen et al. 
19871 iTerasranta etanil992l 1 19981 12004 l2005t iNieppola et al" 
2007, and some unpublished dat a). IR-data from IRA S or 
2MASS catalogues, or from ISO (iPadovani et all l2006t) were 
also available for most sources in the aforementioned databases. 
In addition to X-ray datapoints fo und in the dat a bases , we used 
1 keV monochromatic fluxes from lDonato et al.l (|2001|) . 

The SEDs were plotted in log v - log (vFv) parameter 
space and are presented in Fig.|5](published electronically). The 
synchrotron components of the SEDs were fitted by a simple 
parabolic function 

log(yF,) =A(logv)2 + B(logy) + C, (1) 

where A, B, and C are constants. Thus, the synchrotron peak 
frequency is given by 

logv„ = — . (2) 

With knowledge of Vp, we can calculate the peak flux density, 
VpFyp, using Eq.[T] The peak frequencies obtained in this work 
are listed in Col. (7) of Table[T] 

In Fig. |5] datapoints included in the fitting are marked by 
filled circles and other datapoints are marked by open rectan- 
gles. For the sake of consistency, we excluded datapoints below 
1 GHz from all fits. Below that frequency, the SEDs often fea- 
ture an additional flattening originating from an extended com- 
ponent; we therefore applied a cutoff to avoid any contamination 
in the jet synchrotron fit. The synchrotron peak frequency could 
be determined for all but two sources: 1019-1-222 and 1845H-797. 
These objects have peculiar spectral shapes which could not be 
fitted accurately using a parabolic function. 

Our sample consists of radio-bright AGN, which typically 
have a relatively low synchrotron peak frequency. Thus, in most 
SED fits, we did not include the X-ray datapoints in the syn- 
chrotron component. Our decision was made based on a careful 
visual inspection and the X-ray spectral index when possible. 
In some powerful BLOs, such as 0219-1-428 and 1652-H398, the 
synchrotron component extends to the X-ray domain. In the case 
of 0735+178, the X-ray data are also included in the synchrotron 
component in the interest of constraining the fit and providing a 
more accur ate estimate o f Vp. Th e X-ray indices of the SEDs are 
taken from 'D onato et al.l (1200 lb . For clarity, they are plotted in 
the energy range 1-10 keV and are n ormalized with respe ct to 
the 1 keV monochromatic flux listed in lDonato et al.l (1200 ll) . 

' http://cats.sao.ru 

- http://nedwww.ipac.caltech.edu/ 
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Fig. 1. The variability Doppler factors Dyar plotted against the 
synchrotron peak frequency. The vj, -values are D-corrected. 



There was one case, namely 0430+052, where the syn- 
chrotron peak was exceptionally well sampled, but the parabolic 
function did not fit the data very well. To avoid using an obvi- 
ously flawed Vp, we assigned a new value using the true peak 
defined by the highest point of the synchrotron bump. Thus, we 
obtained log Vp - 13.95, instead of log Vp = 15.94 provided by 
the fit. The same approach was considered for 03 16+41 3, but the 
two methods yielded almost identical results. In the analyses, we 



used the fit value log Vp - 14.40. 



3.2. Doppler factors 

To calculate the de-boosted, intrinsic values of 



and 



VpLp 



we need to calculate the Doppler boosting fac- 
tors D. These have been calculated in the past using 



several methods (e.g. iGhiseUini et al.l 19931: Gtiiiosa & Dal' 
1996t iG uerra & Dalv Il997t iLahteenmaki & Valtaoial 1199' 
Jorstad et al.ii2005tiWu et al.l2007 l). In this work, we use the D^ar 
-values of Hovatta et al. 2008 (in preparation), which have been 
determined fr om total flux density variability data in the same 
manner as in iLahteemnaki & Valta oia (1999)- The flux curves 
from long-term monitoring at 22 and 37 GHz at Metsahovi 
Radio Observatory are decomposed into exponential flares. The 
exponential fit provides the necessary parameters for the calcu- 
lation of the variability brightness temperature ,,0^; comparing 
Tb,var with the eq uipartition value T^a then p rovide s the amount 
of boosting. Both ILahteenmaki & Valtaoial (1 19991) and Hovatta 
et al. (2008) used the equipartition value of 5 x 10'° K proposed 
by^eadhead (1994). 

We were able to derive D^ar for 89 sources in our sample. 
The remainder of the sample had to be excluded from the study 
of intrinsic properties. The values of D^ar and a detailed expla- 
nation of their determination will be published in a future paper 
(Hovatta et al. 2008). 



the redshift-corrected VpLp versus Vp sequence, which has been 
considered in previous papers. We therefore checked how mea- 
surements of D„,r from Hovatta et al. 2008 correlate with the 
D-corrected Vp. The Doppler correction of peak frequency was 
performed using equation 



1 +z 

''Id 



(3) 



where the primed symbol represents the parameter value in the 
source rest frame. The result is presented in Fig. [1] The nega- 
tive correlation is obvious, and the Spearman rank correlation 
test provides a correlation coefficient p - -0.698 and a proba- 
bility P < 0.001 of no correlation. The correlation remains sim- 
ilar if we use values of Vp uncorrected for the Doppler factor. 
The sources with their synchrotron peak at low energies are sig- 
nificantly more boosted than high-v^ sources. A similar nega- 
tive correlation (p - -0.6 54 and P - 0. 004) can be obtained us- 
ing the Doppler factors of I Jorstad et al.l (l2005h . determined from 
variab ility data. We measure d this also using the Doppler fac- 
tors of lGhiselUni etal] (1 19931) . calculated from the difference be- 
tween the observed X-ray flux and one predicted by synchrotron 
self-Compton theory, based on the observed VLBI data. In this 
case, the correlation is very strong (p — -0.517 and P < 0.001). 
Therefore, we conclude that dependence is similar irrespective 
of the method used to calculate the D -factors, and the negative 
correlation between D^.ar and Vp is robust. 

Figure [1] also depicts the broadness of the range of D^,ar, 
stretching from close to to 30. At low Vp in particular, it is cru- 
cial to have an accurate measurement of the source Do ppler fac- 
tor to obtain credible results of the intrinsic properties. I Wu et aTl 
(l2007b determined Doppler factors derived from the relation be- 
tween the observed 5 GHz core luminosity and a theoretical one, 
calculated from the total radio power at 408 MHz. Using log Vp 
-values from Paper I, they also plotted the dependence between 
D and Vp. Their Fig. 13 bears a striking resemblance to Fig.[T]of 
this work. 



4.2. VpLp vs. Vp 



-correlation 



We examined the VpLp versus Vp -correlation using first redshift- 
corrected values and then Doppler-corrected, intrinsic values. 
We observe the results in Fig.|2l In the top panel, we display the 
observational, z-corrected correlation, and in the bottom panel 
the intrinsic, D-corrected correlation. The frequency redshift- 
correction is completed using the equation 



(4) 



and the redshift-corrected luminosity . is calculated with the 
equation 



And] 
r - Lp 



(5) 



where di denotes luminosity distance and Fyp is determined 
from the SED. The Doppler-corrections are performed using 
Eq.|3]and the equation 



4. Results 

4. 1 . The dependence of Dyar and Vp 

For studies of the blazar sequence, the crucial question is 
whether D correlates with Vp and thus changes the shape of 



L' = L„., 



Ill 



assuming F ccv 



(6) 



In all equations, the subscript z represents redshift-corrected val- 
ues and primed values are both redshift and Doppler-corrected. 
When using Eq.|6] we assumed a spectral index a - \ according 
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Fig. 2. The peal: luminosity of the synchrotron component plot- 
ted versus the synchrotron peak frequency. In the top panel, 
the values are z-corrected, and in the bottom panel Doppler- 
corrected. AGN classes are marked with different symbols (see 
legend). 



to the definition of the synchrotron peak of the SED. All equa- 
tions co ncerning redshift- an d Doppler-corrections can be found 
in e.g. lKembhavi & Narlikadd 19991) . 

The blazar sequence is represented by a negative correlation 
between the synchrotron peak frequency Vp and peak luminosity 
VpLp. When our observational, redshift-corrected correlation is 
considered, we find that a significant correlation is also measured 
for our sample. The Spearman rank correlation test measures p 
= -0.248 and P - 0.002 for the entire sample. The Spearman test 
results for different AGN subclasses can be seen in Table|2] only 
the correlation of LPQs is found to be significant. Considering 
the entire sample, the trends visible in the top panel of Fig.|2]are 
similar to those of the blazar sequence scenario found in several 
previous papers. 

The full sample of our study includes many sources that are 
not traditionally classified as blazars, although the blazar classi- 
fication is not clearly defined. For consistency, we identified the 
bona fide blazars in our sample and measured the correlation for 
their data alone. We classified a source as a blazar if it is included 
in one of the following: i) the total blazar sample of Fossati et al. 



1998 ), ii) the extragalactic radio sources list of Wall & Peacock 
19851) and has ajjs ghz < 0.5 (S oc v"), iii) the DXRBS 
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Fig. 3. The Doppler-corrected synchrotron peak luminosity vs. 
peak frequency correlation for bona fide blazars. 



blazar sample of lPerlman et alj (Il998l) : or is classified as a BLO 
according to Table[Tl There were 65 blazars that fulfilled at least 
one of these criteria, and they are marked in Col. (9) of Table[T] 
From Table|2] we observe that the blazar sequence can be found 
in our data when observational quantities are used, because there 
is a significant anticorrelation between the z-corrected VpLp and 
Vp for blazars alone. We list results also for the complete 1.2 
Jy AGN subsample and separately for blazars belonging to the 
1.2 Jy sample. For both subsets of data separately, the negative 
correlation is clearly significant. 

We note that the observational anticorrelation detected here 
is not in disagreement with the conclusion of Paper I, which re- 
ported a lack of correlation between VpLp and Vp for BLOs. The 
range of Vp was very extended in Fig. 6 of that paper, whereas in 
Fig.|2]here the highest value of log Vp is 17.2. When the range 
of log Vp < 17 alone is considered, the two figures become con- 
sistent. We note also that Paper I considered BLOs alone, and no 
redshift correction was applied to their data. 

In the bottom panel of Fig. |2] we observe the effect of ap- 
plying the Doppler-correction to the correlation. It is evident 
that the correlation is altered dramatically. For the entire sam- 
ple, the results of the Spearman test are p = 0.353 and P < 
0.001, which implies that the correlation is positive. In Table|2] 
it is evident that, for BLOs, there is a remarkably clear posi- 
tive correlation between VpLp and Vn. The correlations can be 
interpreted by comparison with Fig. [T] The sources of low Vp are 
more boosted, and their luminosity is lowered more significantly 
by the D-correction. For the blazar sample alone (Fig.|3]l, we also 
find a significant positive correlation (p = 0.366 and P = 0.003). 
Thus, we conclude that the blazar sequence, when defined as an 
anticorrelation of VpLp and Vp, is an artefact of variable Doppler 
boosting across the peak frequency range. 

Interestingly, the positive correlation between the D- 
corrected VpLp and Vp for BL Lac objects can explain the pecu- 
liar U-shape in the VpLp versus Vp -correlation of Paper I (their 
Fig. 6). From Fig. [T] we deduce that for log Vp > 15, D^ar is 
probably s mall, and the_D-correction has little effect (see also 
Fig. 13 of IWu et all (IIoOTI) ). With this in mind, we considered 
datapoints with both log > 15 and a measured redshift from 
Paper I, and applied a z-correction to the Vp-values. We then 
plotted these corrected data together with the D-corrected dat- 
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Table 2. Spearman rank correlation test results for the whole sample and AGN subclasses for the log VpLp vs. log Vp -correlation. 
On the left side is the z-corrected correlation, and on the right side the D-corrected one. Galaxies (GAL) have been omitted due to 
small sample size. 



z-corrected Z)-corrected 



Class 


P 


P 


Class 


P 


P 


All 


-0.248 


0.002 


All 


0.353 


< 0.001 


BLO 


-0.220 


0.126 


BLO 


0.642 


< 0.001 


HPQ 


-0.049 


0.390 


HPQ 


0.354 


0.041 


LPQ 


-0.362 


0.019 


LPQ 


0.236 


0.128 


QSO 


0.110 


0.296 


QSO 


0.550 


0.063 


Blazars 


-0.329 


0.004 


Blazars 


0.366 


0.003 


1.2 Jy AGN sample 


-0.215 


0.013 


1.2 Jy AGN sample 


0.340 


0.002 


1.2 Jy blazars 


-0.287 


0.020 


1.2 Jy blazars 


0.353 


0.009 



apoints for BLOs derived in this paper, producing the top panel 
of Fig. m After combining BLO data sets this way, we are able 
to observe an indisputable positive correlation (p = 0.637 and 
P < 0.001). It appears that the "intrinsic BLO sequence" im- 
plies that sources with higher Vp have higher VpLp: this is com- 
pletely the opposite conclusion to that of the original scenario. 
A similar trend was observed in two BLOs, 1101-H384 (Mrk 
421) and 165 2+398 (Mrk 501), during flares dPian et al.|[T998t 
iTakahashi et al. 2000). Their synchrotron peaks were shifted to 
higher frequencies by the luminosity increase triggered by the 
flare, in contrast to typical observations of the blazar sequence. 
For the case of 1652-1-398 in particular the shift was very pro- 
nounced. The positive correlation between peak frequency and 
luminosity is also in agreement with the number counts of low- 
energy BLOs (LBLs) and high-energy BLOs (HBLs), which in- 
di cate that HBLs a re far less numerous than LBLs. As discussed 
in lPadovani et alj ( l2007h . this is difficult to understand if HBLs 
are intrinsically less luminous than LBLs. 

As a cross-check, we plotted the BLO correlation using in 
addition the -correction, that is the exponent 2 -i- ff in Eq. |6] 
(bottom panel of Fig.|4]l. Both the content of the figure and the 
correlation coefficients of the combined BLO sequence changed 
only marginally. We also checked the D -corrected results of 
Table |2] using the -correction. The significance of the posi- 
tive correlations for the entire sample and the low-Vp BLOs sep- 
arately persists, while for HPQs, LPQs, and QSOs, P > 0.05. 
The main result of this paper, the disappearance of the blazar 
sequence with Doppler-correction, is therefore unchanged. 



5. Discussion 

As described in Paper I, the parabolic SED fit can overestimate 
the synchrotron peak frequency. This occurs in particular when 
the synchrotron peak occurs at or close to X-ray frequencies, 
when no datapoints constrain the parabolic fit beyond the X- 
ray domain. Here we have considered radio-bright quasars and 
blazars, whose peak frequencies are low; in this case, a parabolic 
fit is appropriate and negligible overestimation of the peak fre- 
quency occurs. 

Th e intrinsic peak frequencies have been studied by'Wu et al.' 
(I2007h . who also found a significant negative correlation be- 
tween D and v'^. This independently corroborates the small 
Doppler factors of high-energy BLOs, and provides support to 
the concept of a BLO luminosity sequence presented in Fig. H) 
After all, the log Vp > 15 -datapoints in their Fig. 13, for which 
D < 5 approximately, are exactly the same ones, ta ken from 
Paper I, than those in the high-v^ end of our Fig. H |m et alj 
(l2007h also measured a negative correlation between 408 MHz 
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Fig. 4. Synchrotron peak luminosity versus peak frequency cor- 
relation for BLOs. The data are combined from this work (open 
circles) and Paper I (filled circles). The former are D-corrected, 
the latter are only z-corrected (i.e., we assume that they are not 
Doppler-boosted). In the top panel, the datapoints represented 
by open circles are -corrected, in the bottom panel they are 
-corrected. See text for details. 

intrinsic luminosity and Doppler-corrected Vp. In our view this 
does not, however, correspond to a blazar sequence, which re- 
quires an anticorrelation between Vp and the luminosity at the 
synchrotron peak, VpLp, rather than for any particular frequency 
band. The correlation between radio luminosity and Vp can be 
generated by the shifting of the synchrotron parabola along the 
log V -axis (from low-energy synchrotron source to high-energy 
synchrotron source), regardless of whether a genuine blazar se- 
quence exists. 
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A strong positive correlation between VpLp and Vp for BLOs 
(Fig. was unexpected. For these sources, the cooHng effects 
may decrease and intrinsic luminosity increase as Vp increases. 
However, this phenomenon may not be restricted to BLOs: simi- 
lar correlations may exist for other AGN classes. Unfortunately, 
to our knowledge there are no sufficiently large SED data sets 
that would allow us to investigate their respective luminosity se- 
quences as extensively as that of BLOs. Looking at Table |2l we 
find that HPQs have a marginally significant positive intrinsic 
correlation, and QSOs are close to the limit. With more data- 
points for high values of Vp the correlation might also be sig- 
nificant for QSOs. LPQs, however, appear to differ from other 
classes, exhibiting no intrinsic correlation at all. 

Ano ther open que s tion is the Compton dominance of low-Vp 
sources. iFossati et all ( Il998l) and iGhisellini et all ( Il998l) found 
that as Vp decreased, the inverse Compton component became 
more luminous compared with its synchrotron counterpart in the 
SED. The accurate, observational determination of the shape and 
size of the Compton component for a large sample of objects 
would require a large amount of data from X-ray to TeV ener- 
gies. Such a database is unfortunately not yet available for BLOs 
nor other AGN. 

6. Conclusions 

We have studied the correlation between the synchrotron peak 
frequency, Vp, and the synchrotron peak luminosity, VpLp, and 
the way that this changes when a Doppler-correction is applied 
to both quantities. We have used an extensive sample of 135 
radio-bright AGN, a large part of which forms a complete flux- 
limited 1.2 Jy northern AGN sample. Our conclusions are as fol- 
lows: 

1. There is a clear negative correlation between the Doppler 
factor, D, and Vp, independent of the method for calculating 
D. This correlation has been tested using Doppler factors cal- 
culated from variability data, core versus extended flux, and 
synchrotron self-Compton flux. Typically, for sources with 
log Vp > 15, D is small, while low-Vp sources may have val- 
ues of D as high as 30 and thus be significantly boosted by 
comparison. 

2. The negative correlation between VpLp and Vp at low Vp, 
known as the blazar sequence, is not present when the intrin- 
sic, Doppler-corrected values are used. Instead, the correla- 
tion becomes positive. This is true for the entire AGN sample 
and a subsample of bona fide blazars. The blazar sequence 
is an artefact of variable Doppler boosting across the peak 
frequency range. 

3. For BLOs in particular, there appears to be a strong intrinsic 
positive correlation between VpLp and Vp. This correlation 
is further strengthened when the BLO sample in this work 
is combined with the BLO sample of Nieppola et al. ( 200^. 
The positive correlation is also in accordance with number 
counts of LBLs and HBLs. 
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Table 1. The sample and the observational synchrotron peak frequencies obtained in this work. The values of Vp listed here have been determined 
directly from the SEDs and are not corrected for redshift nor for Doppler boosting. 



Source 


Alias 


R.A.(J2000) 


Dec(J2000) 


1.2 JyAGN sample 


Class 


logvp 


z 


Blazar sample 


0003-066 


NRA0 5 


00:06:13.90 


-06:23:35.34 




BLO 


12.80 


0.347 


* 


0007+106 


PG 0007+106 


00:10:31.01 


+ 10:58:29.50 


* 


GAL 


13.34 


0.089 




0014+813 




00:17:08.47 


+81:35:08.10 


* 


QSO 


13.91 


3.366 




0016+731 




00:19:45.79 


+73:27:30.02 


* 


LPQ 


12.82 


1.781 




0048-097 


PKS 0048-097 


00:50:41.32 


-09:29:05.21 




BLO 


13.28 


0.300 


* 


0059+581 




01:02:45.76 


+58:24:11.14 


* 


QSO 


12.49 


0.644 




0106+013 


OC012 


01:08:38.77 


+01:35:00.30 




HPQ 


12.73 


2.107 




0109+224 


S2 0109+22 


01:12:05.80 


+22:44:39.00 


* 


BLO 


13.81 




* 


0133+476 


DA 55 


01:36:58.59 


+47:51:29.10 


* 


HPQ 


12.77 


0.860 


* 


0149+218 




01:52:18.06 


+22:07:07.70 


* 


LPQ 


13.04 


1.320 




0202+149 


4C 15.05 


02:04:50.41 


+ 15:14:11.00 


* 


HPQ 


12.28 


0.405 


* 


0212+735 




02:17:30.81 


+73:49:32.60 


* 


HPQ 


12.34 


2.367 


* 


0215+015 




02:17:48.95 


+01:44:49.70 


* 


HPQ 


13.49 


1.715 




0218+357 


S4 0218+357 


02:21:05.40 


+35:56:15.00 


* 


QSO 


12.75 


0.936 




0219+428 


3C 66A 


02:22:39.60 


+43:02:08.00 


* 


BLO 


14.96 


0.444 


* 


0224+671 




02:28:50.05 


+67:21:03.00 


* 


QSO 


12.66 


0.523 




0229+131 




02:31:45.89 


+ 13:22:54.72 


* 


LPQ 


13.01 


2.059 




0234+285 


4C 28.07 


02:37:52.41 


+28:48:09.00 


* 


HPQ 


12.81 


1.207 




0235+164 


AO 0235+164 


02:38:38.80 


+16:36:59.00 


* 


BLO 


13.31 


0.940 


* 


0248+430 




02:51:34.54 


+43:15:15.83 


* 


LPQ 


13.45 


1.311 




0300+470 


4C 47.08 


03:03:35.20 


+47:16:17.00 


* 


BLO 


14.34 


0.475 


* 


0306+102 


PKS 0306+102 


03:09:03.60 


+ 10:29:16.00 


* 


BLO 


12.79 


0.863 


* 


0316+413 


3C 84 


03:19:48.16 


+41:30:42.10 


* 


GAL 


14.40 


0.017 


* 


0326+349 




03:29:15.35 


+35:10:08.10 


* 




12.99 


0.500 




0333+321 


NRAO 140 


03:36:30.11 


+32:18:29.30 


* 


LPQ 


13.28 


1.259 




0336-019 


CTA 026 


03:39:30.94 


-01:46:36.00 




HPQ 


13.09 


0.852 




0355+508 


NRAO 150 


03:59:29.75 


+50:57:50.20 


* 


QSO 


12.53 


1.510 




0415+379 


3C 111 


04:18:21.28 


+38:01:35.80 


* 


GAL 


13.64 


0.049 




0420-014 


OA 129 


04:23:15.80 


-01:20:33.10 




HPQ 


12.87 


0.915 


* 


0422+004 


OF 038 


04:24:46.80 


+00:36:07.00 




BLO 


14.83 


0.310 




0430+052 


3C 120 


04:33:11.10 


+05:21:15.60 


* 


GAL 


13.95' 


0.033 


* 


0440-003 


NRAO 190 


04:42:38.66 


-00:17:43.00 




HPQ 


13.50 


0.844 


* 


0446+112 


PKS 0446+112 


04:49:07.67 


+ 11:21:28.00 


* 


GAL 


12.58 


1.207 




0454+039 




04:56:47.10 


+04:00:53.00 




LPQ 


14.11 


1.343 




0458-020 


PKS 0458-020 


05:01:12.81 


-01:59:14.00 




HPQ 


12.62 


2.286 




0507+179 




05:10:02.37 


+ 18:00:41.58 


* 


QSO 


12.71 


0.416 




0528+134 


PKS 0528+134 


05:30:56.42 


+13:31:55.15 


* 


LPQ 


12.45 


2.070 


* 


0552+398 


DA 193 


05:55:30.81 


+39:48:49.17 


* 


LPQ 


12.12 


2.365 




0605-085 


PKS 0605-085 


06:07:59.70 


-08:34:50.00 




HPQ 


13.40 


0.872 


* 


0642+449 


OH 471 


06:46:32.03 


+44:51:16.59 


* 


LPQ 


12.39 


3.406 




0716+714 


S5 0716+714 


07:21:53.30 


+71:20:36.00 


* 


BLO 


13.87 


0.300 


* 


0723-008 


PKS 0723-008 


07:25:50.64 


+00:54:56.54 




GAL 


14.39 


0.127 




0735+178 


PKS 0735+17 


07:38:07.40 


+ 17:42:19.00 


* 


BLO 


14.07 


0.424 


* 


0736+017 




07:39:18.03 


+01:37:04.60 


* 


HPQ 


13.42 


0.191 


* 


0738+313 




07:41:10.70 


+31:12:00.23 


* 


GAL 


13.54 


0.631 




0748+126 




07:50:52.05 


+ 12:31:04.83 


* 


LPQ 


12.76 


0.889 




0754+100 


PKS 0754+100 


07:57:06.64 


+09:56:34.90 


* 


BLO 


13.70 


0.266 


* 


0804+499 




08:08:39.67 


+49:50:36.50 


* 


HPQ 


12.75 


1.430 




0814+425 


OJ425 


08:18:16.00 


+42:22:45.40 


* 


BLO 


13.26 


0.258 


* 


0823+033 


PKS 0823+033 


08:25:50.30 


+03:09:24.00 


* 


HPQ 


12.73 


0.506 


* 


0827+243 


OJ 248 


08:30:52.08 


+24:10:60.00 


* 


LPQ 


12.94 


0.941 




0829+046 


PKS 0829+046 


08:31:48.88 


+04:29:39.09 


* 


BLO 


13.92 


0.180 


* 


0836+710 


4C 71.07 


08:41:24.37 


+70:53:42.20 




LPQ 


13.94 


2.170 




0846+513 




08:49:57.98 


+51:08:29.00 




QSO 


13.24 


1.860 




0850-1213 




08:50:09.60 


-12:13:34.00 




QSO 


12.58 


0.566 




0851+202 


OJ287 


08:54:48.80 


+20:06:30.00 


* 


BLO 


13.64 


0.306 


* 


0906+430 


3C216 


09:09:33.50 


+42:53:46.08 


* 


HPQ 


12.99 


0.670 


* 


0917+449 




09:20:58.46 


+44:41:53.99 


* 


QSO 


12.85 


2.180 




0923+392 


4C 39.25 


09:27:03.01 


+39:02:20.90 


* 


LPQ 


11.99 


0.699 


* 


0945+408 


4C 40.24 


09:48:55.34 


+40:39:44.60 


* 


LPQ 


12.87 


1.252 
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Source 


Alias 


R.A.(J2000) 


Dec(J2000) 


1 .2 Jy AGN sample 


Class 


log Vp 


z 


Blazar sample 


0953+254 




09:56:49.88 


+25:15:16.05 




LPQ 


13.03 


0.712 




0954+556 


S4 0954+556 


09:57:38.24 


+55:22:58.00 


* 


HPQ 


13.77 


0.901 


* 


0954+658 


S4 0954+65 


09:58:47.20 


+65:33:54.00 


* 


BLO 


14.16 


0.367 


* 


1019+222 




10:21:54.59 


+21:59:31.20 


* 


GAL 




1.617 




1036+054 




10:38:46.78 


+05:12:29.09 






12.57 


0.473 




1040+244 


TEX 1040+244 


10:43:09.00 


+24:08:35.00 


* 


BLO 


12.82 


0.560 




1039+811 




10:44:23.06 


+80:54:39.44 


* 


LPQ 


13.04 


1.260 




1049+215 




10:51:48.79 


+21:19:52.31 




LPQ 


12.84 


1.300 




1055+018 


OL093 


10:58:29.61 


+01:33:58.80 




HPQ 


12.68 


0.888 




1150+812 




11:53:12.50 


+80:58:29.15 




QSO 


12.75 


1.250 




1150+497 




11:53:24.47 


+49:31:08.83 


* 


HPQ 


13.29 


0.334 




1156+295 


4C 29.45 


11:59:31.83 


+29:14:44.00 




HPQ 


13.44 


0.729 




1219+285 


ON 231 


12:21:31.70 


+28:13:58.00 




BLO 


15.10 


0.102 


* 


1219+044 




12:22:22.55 


+04:13:15.78 




QSO 


13.07 


0.965 




1222+216 


PKS 1222+216 


12:24:54.51 


+21:22:47.00 


* 


QSO 


13.27 


0.435 




1226+023 


3C 273 


12:29:06.69 


+02:03:08.60 


* 


LPQ 


13.89 


0.158 


* 


1253-055 


3C 279 


12:56:11.17 


-05:47:21.50 




HPQ 


12.67 


0.538 


* 


1308+326 


AUCVn 


13:10:28.66 


+32:20:43.80 




BLO 


13.49 


0.992 


* 


1324+224 




13:27:00.86 


+22:10:50.16 




QSO 


12.74 


1.400 




1406-076 


PKS 1406-076 


14:08:56.48 


-07:52:27.00 




QSO 


12.69 


1.494 




1413+135 


PKS 1413+135 


14:15:58.80 


+ 13:20:24.00 


* 


BLO 


12.64 


0.247 


* 


1418+546 


OQ 530 


14:19:46.60 


+54:23:14.00 


* 


BLO 


14.03 


0.152 




1502+106 


OR 103 


15:04:24.98 


+ 10:29:39.00 


* 


HPQ 


13.08 


1.833 




1510-089 


PKS 1510-089 


15:12:50.53 


-09:05:59.00 




HPQ 


13.49 


0.361 




1532+016 


PKS 1532+016 


15:34:52.45 


+01:31:04.00 


* 


HPQ 


13.23 


1.435 




1538+149 


4C 14.60 


15:40:46.50 


+ 14:47:45.90 


* 


BLO 


13.84 


0.605 


* 


1546+027 




15:49:29.44 


+02:37:01.20 




HPQ 


12.97 


0.413 




1548+056 




15:50:35.27 


+05:27:10.45 




HPQ 


12.74 


1.422 




1606+106 


4C 10.45 


16:08:46.20 


+10:29:07.80 




LPQ 


12.80 


1.226 




1611+343 


DA 406 


16:13:41.00 


+34:12:48.00 


* 


LPQ 


12.69 


1.401 


* 


1633+382 


4C 38.41 


16:35:15.49 


+38:08:04.50 




LPQ 


12.89 


1.814 




1637+574 


OS 562 


16:38:13.46 


+57:20:24.00 




LPQ 


13.17 


0.745 




1642+690 




16:42:07.85 


+68:56:39.76 




HPQ 


12.54 


0.751 




1641+399 


3C 345 


16:42:58.81 


+39:48:37.00 




HPQ 


12.71 


0.595 




1652+398 


MRK501 


16:53:52.20 


+39:45:36.00 


* 


BLO 


17.19 


0.034 


* 


1725+044 


PKS 1725+044 


17:28:24.95 


+04:27:05.00 


* 


QSO 


13.27 


0.297 


* 


1730-130 


NRAO 530 


17:33:02.71 


-13:04:49.55 




QSO 


12.42 


0.902 




1739+522 


S4 1739+52 


17:40:36.98 


+52:11:43.00 


* 


HPQ 


12.93 


1.379 




1741-038 


PKS 1741-038 


17:43:58.86 


-03:50:04.60 




HPQ 


12.35 


1.054 


* 


1749+096 


PKS 1749+096 


17:51:32.70 


+09:39:01.00 


* 


BLO 


12.78 


0.320 




1758+388 




18:00:24.72 


+38:48:31.10 




QSO 


12.70 


2.092 




1803+784 


S5 1803+784 


18:00:45.40 


+78:28:04.00 


* 


BLO 


13.60 


0.684 




1807+698 


3C 371 


18:06:50.70 


+69:49:28.00 


* 


BLO 


14.46 


0.051 


* 


1823+568 


4C 56.27 


18:24:07.07 


+56:51:01.50 


* 


BLO 


12.93 


0.663 


* 


1828+487 


3C 380 


18:29:31.80 


+48:44:46.62 




LPQ 


13.88 


0.692 




1845+797 


3C 390.3 


18:42:08.99 


+79:46:17.00 




GAL 




0.055 




1901+319 




19:02:55.93 


+31:59:41.70 




QSO 


15.67 


0.635 




1926+61 


S4 1926+61 


19:27:30.40 


+61:17:32.00 


* 


BLO 


13.27 




* 


1928+738 


4C 73.18 


19:27:48.50 


+73:58:01.60 


* 


LPQ 


13.32 


0.303 


* 


1954+513 




19:55:42.74 


+51:31:48.55 




LPQ 


12.90 


1.223 




2007+776 


S5 2007+77 


20:05:31.10 


+77:52:43.00 




BLO 


12.95 


0.342 


* 


2005+403 




20:07:44.94 


+40:29:48.60 




QSO 


13.07 


1.736 




2021+614 


OW637 


20:22:06.68 


+61:36:58.80 


* 


LPQ 


12.67 


0.227 


* 


2022+171 




20:24:56.56 


+ 17:18:13.20 




LPQ 


13.14 


1.050 




OA'2T 1^11 




ZU:jo: J /.U4 


1 ^ 1 . 1 0. 1 O 

+jl:iy:lz.oo 




QoU 


iz.zo 


l.OO / 




2059+034 




21:01:38.83 


+03:41:31.30 


* 


QSO 


13.03 


1.013 




2121+053 




21:23:44.52 


+05:35:22.09 


* 


HPQ 


13.01 


1.941 




2134+004 


OX 057 


21:36:38.59 


+00:41:54.21 


* 


LPQ 


12.32 


1.936 


* 


2136+141 




21:39:01.31 


+ 14:23:36.00 


* 


LPQ 


12.71 


2.427 




2141+175 




21:43:35.55 


+ 17:43:48.50 


* 


LPQ 


13.70 


0.211 




2144+092 




21:47:10.16 


+09:29:46.67 




QSO 


13.04 


1.113 





Table 1. continued. 
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Source 


Alias 


R.A.(J2000) 


Dec(J2000) 


1.2 Jy AGN sample 


Class 


log 


z 


Blazar sample 


2145+067 




21:48:05.46 


+06:57:38.60 




LPQ 


12.73 


0.990 




2200+420 


BL LAC 


22:02:43.30 


+42:16:39.00 




BLO 


14.14 


0.070 




2201+315 


4C 31.63 


22:03:14.98 


+31:45:38.30 




LPQ 


13.63 


0.298 




2201+171 




22:03:26.89 


+ 17:25:48.20 


H: 


QSO 


12.90 


1.075 




2209+236 


PKS 2209+236 


22:12:05.97 


+23:55:40.00 




QSO 


12.62 


1.125 




2223-052 


3C446 


22:25:47.30 


-04:57:01.39 




BLO 


12.74 


1.404 




2227-088 




22:29:40.09 


-08:32:54.50 




HPQ 


12.99 


1.562 








ZZ.jZ.JV.^- L 


+i i.4j.ju.yu 








1 M'iT 
i .Uj / 




2234+282 




22:36:22.47 


+28:28:57.41 


* 


HPQ 


12.74 


0.795 




2251+158 


3C 454.3 


22:53:57.75 


+ 16:08:53.60 


* 


HPQ 


12.83 


0.859 


* 


2254+074 


PKS 2254+074 


22:57:17.30 


+07:43:12.27 




BLO 


14.07 


0.190 


* 


2344+092 




23:46:36.84 


+09:30:45.52 


* 


LPQ 


14.25 


0.677 


* 


2351+456 


4C 45.51 


23:54:21.68 


+45:53:04.00 


* 


QSO 


12.44 


1.986 




2353+816 




23:56:22.79 


+81:52:52.26 


* 


QSO 


12.74 


1.344 





' = log Vp determined from the highest flux density measurement of the synchrotron component of the SED (see Sect. l3.lt . 
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Fig. 5. The spectral energy distributions determined in this work. All quantities are observational and uncorrected for redshift and 
Doppler boosting. The dashed hne represents the paraboUc fit to the synchrotron component. Datapoints marked with filled circles 
are included in the fit, the empty squares are not. Sources 1019-1-222 and 1845-1-797 could not be fitted with a parabolic function. 
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Fig. 5. continued. 
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Fig. 5. continued. 
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Fig. 5. continued. 
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Fig. 5. continued. 
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Fig. 5. continued. 
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Fig. 5. continued. 
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Fig. 5. continued. 



